Análise Matemática IV 

Problemas para as Aulas Práticas 



Semana 10 



1. Determine a solução da equação linear: 

y^ - 2y^ + y'-2 = b(t) 

que verifica as condições iniciais 

y(0) = y'(0) = , y (2) (0) = l 

quando: 

(i) 6(í) = 0, Vi G R. 

(ii) b(t) =t, Vi G R. 

(iii) ò(í) = e*, Vi G R. 

Resolução: 

A solução da equação 

y< 3 > - 2y^ + y' = b(t) + 2 

é 

y(t) = y h (t) + y p (t) 

em que yh(t) é a solução geral da equação homogénea associada, e y p (t) é uma solução 
particular da equação completa. 

(i) Sendo ò(i) = 0, estamos a resolver a equação 

y^ - 2y^ + y' = 2 

e é óbvio que y p (t) = 2t é uma solução particular da equação. Para determinar a solução 
geral da equação homogénea 

y^ - 2y( 2 ) + y> = 

note-se que a equação característica é 

A 3 - 2A 2 + A = ^ A(A - l) 2 = 

Então 

yi(t) = e ot , y 2 (t) = e t , y 2 (t) = té 

são 3 soluções linearmente independentes da equação homogénea, pelo que 

y h (t) = a + be 1 + cté 
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Finalmente, a solução da equação pedida é 

y(t) = y h (t) + y p (t) = a + bê + cte 1 + 2í 
Dado que y(0) = y'(0) = e y"(0) = 1, concluimos que 



' a + 6 = 



a = 5 



<6 + c + 2 = ^><6 = -5 



6 + 2c = 1 



c = 3 



pelo que 

(ii) A solução da equação 



y(t) = 5 - 5e* + Síe* + 2í 
y (3) _ 2y (2) + y < = t + 2 



(1) 



será 



y(í) = y h (t) + y p (t) 

em que yh(t) é a solução geral da equação homogénea associada, e y p (t) é uma solução 
particular da equação completa. Como calculado em (i) 

y h (t) = a + be 1 + cte 1 

Para calcular y p , note-se que b(t) = t + 2 = (í + 2)e OÍ é solução da equação diferencial 

y" = D 2 y = 

Temos então, aplicando P a {D) = D 2 à equação (1) 

D 2 (D 3 - D 2 + D)y = D 2 (t + 2) O L> 3 (L> - l) 2 y = 

Esta equação admite como solução geral (atenda a que as soluções do polinómio carac- 
terístico são 1 de multiplicidade 2 e de multiplicidade 3) 

y(t) = a + bt + et 2 + de t + fte t 

e conclui-se que bt + et 2 é a candidata a solução particular (dado que a + de t + /te* é y^). 
Resta agora calcular as constantes b e c de modo a que y p (í) = bt + et 2 seja solução de 
(1). Assim 



2/ p (í) = 6í + cí 2 , y' p (t) = b + 2ct , y ;'(t) = 2c , y' p "{t) = 



pelo que 



y£ 3) - 2y p 2 > + y ; = í + 2 ^ 2ct + b - Ac = t + 2 O c = - e 6 = 4 



Tem-se então que y p (í) = At + y, e a solução geral de (1) é 



y(t) = a + be* + cte 1 + At + 
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Dado que y(0) = y'(0) = e y"(0) = 1, concluímos que 



( a+b=0 



a = 



<6 + c + 4 = o- < b = 



b + 2c+ 1 = 1 



c = 4 



pelo que 



y(t) = 8 - 8e* + 4íe* + 4í + 



(iii) A solução da equação 



,( 3 ) 



2y 



( 2 ) + y' = e * + 2 



(2) 



será 



í/(í) = j/ h (í) + y p (í) 



em que yh(t) é a solução geral da equação homogénea associada, e y p (t) é uma solução 
particular da equação completa. Como calculado em (i) 

y h (t) = a + be 1 + ctê 

Para calcular y p , começamos por verificar que como consequência da linearidade da deri- 
vada 

y p (t) = 2t + z p (t) 

em que z p (t) é solução particular da equação 



y 



(3) 



2y 



(2) 



y = e 



(3) 



Sendo 6(í) = é solução da equação diferencial (por estar associada à raiz do polinómio 
característico 1 de multiplicidade 1) 

y' - y = {D - \)y = 

teremos, aplicando Pa{D) = D — 1 à equação (3) 

(D - 1)D(D - ify = (D - l)é O D(D - ify = 

Esta equação admite como solução geral (atenda a que as soluções do polinómio carac- 
terístico são de multiplicidade 1 e 1 - de multiplicidade 3) 

y(t) = a + be l + cte 1 + dfe 1 

e conclui-se que dí 2 e í é a candidata a solução particular (dado que a + be l + cte 1 é y^). 
Resta agora calcular a constante d de modo a que z p (t) = dt 2 e t seja solução de (3). Assim 

z p (t) = dt 2 é , z' p (t) = d(2t + í 2 )e* , z' p '(t) = d(2 + At + t 2 )é , z' p "(t) = d(6 + 6í + í 2 )e* 

pelo que 



4 3) - 24 2) + z' = e* O 2d = 1 O d= - 
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Tem-se então que z p (t) = *-£-, e a solução geral de (2) é 



y(t) = a + òe* + cté + 2í + 



tV 



Dado que y(0) = y'(0) = e y"(0) = 1, concluímos 



( a+b=0 f o=4 



<6 + c + 2 = ^ ^ ò = -4 



6 + 2c+ 1 = 1 



c = 2 



pelo que 



y (í) = 4 - Ae t + 2te t + 1t 



íV 



2. Determine a solução da equação diferencial 

y"-4y' + 3y=(l + e- x )- 1 

que verifica as condições iniciais y(0) = y'(0) = 1 
Resolução: 

A solução da equação é da forma 

y(t) = yh(x) + y p (x) 

sendo y^ a solução geral da equação homogénea associada e y p uma solução particular da 
equação completa. Comecemos por calcular y^. O polinómio característico associado é 

A 2 -4A + 3 = O A=louA = 3 

Sendo assim, e x e e 3x são duas soluções linearmente independentes da equação homogénea, 
pelo que 

y h (t) = Cl e x + c 2 e 3x 

Para calcular y p , notamos que teremos forçosamente de utilizar a fórmula da variação das 
constantes, dado que 

h(x) = 

não é solução de nenhuma equação diferencial linear de coeficientes constantes. Assim 
sendo 

r r o 

y p (x) = [e x e 3x ] W^s) , ds 



l+e" 



em que W(x) é a matriz Wronskiana associada, isto é 

W(x) -- 



pX páX 

e x 3e 3x 
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Assim 



y p (x) = [ e x e 3x ] 



= [ e x e 3x ] 



pX p<->X 



3e" 



_p à- s p ■- >s 



O 

i 



l+e- 



ds 






ds 



log(l + e~ x ) 
- + e~ x — log(l + e~ x ) 



2 

\ U log(l + e"*) - Ç + e 2 * - e 3 * log(l + O 






log(l + e" 



2 4 2 

— 3s 

onde calculámos a primitiva de -^—j fazendo a substituição e~ s = t. Finalmente 



y(x) = y h (x) + y p (x) = c x e x + c 2 e Ax + 
Dado que y(0) = y'{0) = 1 tem-se que 

3 






log(l + e~ x ) 



Cl + c 2 = 



ci = l-Mog2 



pelo que 



ci + 3c 2 = | + log2 
y(x) = 



p x p^X 



log 



O 



1 + e"* 



c 2 = -± + ±log2 



4 ' 2 



g 3x g 2x g e a 



3. Considere a equação 



y (3) _ 4y (2) + 5y / = o 

(i) Determine a sua solução geral. 

(ii) Determine para que condições iniciais em í = é que os problemas de valor inicial 
correspondentes têm soluções convergentes quando t —> oo. 

(i) polinómio característico associado é 

A 3 - 4A 2 + 5A = A(A 2 -4A + 5) = <=> \ = 0ou\ = 2±i 

pelo que 1, e 2í cosí e e 2í sení são três soluções linearmente independentes da equação. 
Como tal a sua solução geral é 

y(t) = Ci + c 2 e cosí + C3e sent 

(ii) Pretende-se determinar a, (3 e 7 de modo a que o problema de valor inicial 

3,(3) _ 43/(2) + 5 y ' = o 

3/(0) = a , y'(0) = , y"(0) = 7 
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tenha solução limitada. Observe que a função y(t) possui limite quando t —>■ oo sse 
c 2 = c 3 = 0, isto é se é da forma y(t) = C\. Sendo assim, as condições iniciais requeridas 
serão 

3/(0) = a , y'(0) = , y"(0) = 

para «eR. 

4. Para que valores de c G IR é que a equação 

y" - 2q/ + y = o 

admite uma solução periódica, que não seja identicamente nula? 
Resolução: O polinómio característico associado é 

A 2 - 2cA + 1 = «• A = c ± \/c 2 - 1 

Há, então, três hipóteses: 

• c 2 — 1 = 0: Neste caso o polinómio característico tem apenas uma raiz c pelo que a 
solução da equação diferencial é 

y(t) = Cl e ct + c 2 te ct . 

A única solução periódica é a solução nula: c\ = c 2 = 0. 

• c 2 — 1 > 0: Neste caso o polinómio característico tem duas raízes reais distintas c + 
e c~ , pelo que a solução da equação diferencial é 

y(t) = Cl e c+t + c 2e c_í . 

Mais uma vez, a única solução periódica é a solução nula: c\ = c 2 = 0. 

• c 2 — 1 < 0: Neste caso o polinómio característico tem duas raízes complexas conju- 
gadas c ± i\J\ — c 2 , pelo que a solução da equação diferencial é 

y(t) = c ie ct cos(tVl-c 2 ) + c 2 e ct sen(tVl-c 2 ) 

Para que esta solução seja periódica teremos de ter c = 0. Assim, este valor é o 
único valor de c para o qual a equção admite soluções periódicas não triviais. Neste 
caso as soluções são: 

y(i) = Ci cos(í) + c 2 sen(í) 

5. Considere a equação 

2,(4) + 2y( 3 ) + y( 2 ) = t + COS t (4) 

(i) Determine a solução geral da equação homogénea correspondente a (4). 
(ii) Determine uma solução particular de (4). 
(iii) Determine a solução de (4) que verihca a condição inicial 

j/(0) = y / (0) = y( 2 )(0) = y( 3 )(0) = 



Análise Matemática IV 7 

Resolução: 

(i) polinómio característico associado é 

A 4 + 2A 3 + A 2 = & A 2 (A + l) 2 = 

pelo que 1, t, e _í e te~ t são soluções independentes da equação homogénea, e como tal 

Vh{t) = ci + c 2 t + c 3 e~* + c^te' 1 

é a solução pedida. 

(ii) Como usualmente a solução de (4) é dada por 

y(t) = y h (t) + y p (t) = c x + c 2 t + c 3 e~ t + c^te' 1 + y p (t) 

em que y p é uma solução particular de (4). Para utilizar o método dos coeficientes 
indeterminados, podemos considerar 

y P {t) = z x {t) + z 2 {t) 

em que z\ é solução da equação diferencial 

yW + 2y^ + y^ = t (5) 

e, z 2 é solução da equação diferencial 

y (4) + 2y (3) + y {2) = cosi (6) 

Comecemos por calcular z\. Visto t ser solução da equação 

y" = D 2 y = 

aplicando P a (D) = D 2 à equação (5) 

D A (D + l) 2 y = D 2 t = 

Esta equação admite como solução geral (atenda a que as soluções do polinómio carac- 
terístico são de multiplicidade 4 e -1 de multiplicidade 2) 

y(t) = a\ + a 2 t + a 3 t 2 + a^t 3 + a^e~ l + a^te -1 

e conclui-se que a^ + a^t 3 é a candidata a solução particular (dado que ai + a 2 t + az,e~ t + 
a^te -1 é yh). Resta agora calcular as constantes a 3 e a^ de modo a que z\(t) = a 3 t 2 + a^t 3 
seja solução de (5). Assim 

zi(í) = a 3 t 2 +a A t 3 , z[(t) = 2a 3 t+3a A t 2 , z"(t) = 2a 3 +6a 4 t , z'"(t) = 6a 4 , z^°\t) = 
pelo que 

4 4) + 2zf } + z[ 2) = t & 12a 4 + 2a 3 + Ç>a A t = t <& \ az ^~ l 
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Tem-se então que Zi(t) = — t + 1-. Para calcular z 2 , note-se que cosi é solução da equação 

y" + y = {D 2 + l)y = 

aplicando P a (D) = D 2 + 1 à equação (6) 

(D 2 + 1)D 2 (D + l) 2 y= (D 2 + l)cosí = 

As raízes do polinómio característico são (com multiplicidade 2), -1 (com multiplicidade 
2) e ±z (simples). Assim, a equação acima admite como solução geral: 

y(t) = cii + a 2 t + a 3 e _í + a 4 íe _í + a 5 cos t + a 6 sen £ 

e conclui-se que a 5 cosí + a 6 sení é a candidata a solução particular (dado que a^ + 
«2^ + «3e _í + a 4 íe _í é y^,). Resta agora calcular as constantes a*, e ci6 de modo a que 
2:2 (í) = «5 cosi + a^sent seja solução de (6). Assim 

z 2(t) = «5 cosi + as sení , 2 2 (í) = — 05 sení + aç cosi , z^t) = —a^ cos t — a^ sen t 
z' 2 "{t) = a 5 sent — a e cost , z 2 (t) = z 2 (t) 



pelo que 

z y 2 ' + 2z K 2 ' + z K 2 ' = t <=? —2ae cosi + 2a$ sen t = cosi -£4> 
Tem-se então que 2 2 (£) = — \ sen í. Finalmente 



a 5 = 
a 6 = - 



y(í) = s/ fc (í) + zi(t) + z 2 (t) = ci + c 2 í + c 3 e * + c 4 íe * - t 2 
(iii) Dado que y(0) = y'(0) = y"(0) = y"'{0) = concluimos 



sení 



ci + c 3 = 
c 2 — c 3 + C4 — 2 = 
c 3 - 2c 4 - 2 = 
-c 3 + 3c 4 + | = 



<^ < 



ci = — 3 

c 2 = 3 

c 3 = 3 

c 4 = 5 



pelo que 



y(í) = -3 + 3í + 3e _í + -te 



sení 



